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Astronomical observations of elemental and isotopic abundances in a variety of environments 
provide the means to extract the source and production route for elements as a function of 
time since our Galaxy formed about 15 billion years ago. Spallation reactions, where Galactic 
cosmic rays break apart interstellar nuclei, play an important role in the synthesis of the light 
elements - lithium, beryllium, and boron ^. Measurements of the interstellar lithium isotope 
ratio offer a test of the proposed models of lithium synthesis in particular and light element 
synthesis in general Here we report the presence of variations in this ratio in interstellar 
cloud(s) toward stars associated with a site of massive star formation. Our observations 
reveal a sight line (o Persei) with ''Li/^Li ratios of approximately 2, rather than the Solar 
System ratio of about 12.3 While our measured value of approximately 2 provides the 
first clear evidence for undiluted products of cosmic ray spallation in interstellar space and 
confirms a basic tenet of models of light element synthesis the expected enhancement in 
elemental lithium abundance is not seen. 

The observations were acquired with the Harlan J. Smith 2.7 meter telescope and the 
"2dcoude" spectrograph at McDonald Observatory We observed ( Per (Bl lb; V = 2.85; 
V sini = 59 km s~^) in 1996 November and 1997 January and o Per (Bl HI; V = 3.83; v sini 
= 85 km s~^) in 1998 January. High resolution spectra (R ~ 170,000) in two spectral regions 
were obtained, at 4044 A for K I and 6707 A for Li I. Figure 1 displays the spectra, after 
removal of the undulating stellar continuum through the use of a low order polynomial. 

Generally there are multiple interstellar clouds on any given line of sight whose existence 
is revealed by absorption at distinct Doppler velocities. The velocity structure has to be 
well known before an accurate lithium isotope ratio can be deduced from the spectra 
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Therefore, the observed directions should have relatively simple velocity structure (i.e., one 
or two interstellar clouds with detectable Li I lines). In addition, data on another species, 
which has similar properties to lithium and resides in the same portion of the interstellar 
cloud, should be obtained for use as a template of the velocity structure. Our method 
differs slightly from the technique of Lemoine et al. [Ref. 9, 10]. Lemoine et al. used K I 
A7699 as their velocity template. This line of K I is much too strong to use as an effective 
velocity template. Instead, we obtained data on K I A4044 which is of comparable strength 
to the lithium lines (see Fig. I) and therefore a more appropriate template. (The oscillator 
strength for A4044 is about a factor of 50 smaller ^^.) Moreover, high signal to noise, high 
resolution spectra of the Li doublet at 6707 A are required, since the lines are weak and 
the fine structure splitting for ''Li I is comparable to the isotope shift of ^Li I (~ 0.160 A), 
resulting in a blend of the ''Li and ^Li lines as shown in Fig. I. 

The observations were fitted by synthetic profiles where profile parameters were adjusted 
until a minimum was obtained; this method was successfully applied in our determination 
of the interstellar ^^B/^°B ratio i2,i3_ ^pj-^g results from applying the parameters obtained 
independently from the ''Li I and K I lines show agreement in column densities at the 1-a 
level. The results of profile fitting are displayed in Table I, and Table 2 contains the inferred 
isotope ratios. 

Our ''Li/^Li ratio of 10.6 it 2.9 for the cloud toward ^ Per is just consistent with a 
previous determination (5.5ti'i), the difference being in the derived A^(^Li I), and agrees 
with the most recent meteoritic value '' (II. 86-12. 44). As shown by inspection of the line 
profile, the ''Li/^Li ratios for o Per's clouds are highly non-solar: ''Li/^Li = 3.6 it 0.6 and 
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1.7 lb 0.3, for the stronger and weaker cloud, respectively, with the latter being essentially 
the value (~ 2) predicted for Li production by cosmic ray-induced spallation. The most 
recent models "^"^ use up-to-date cross sections for spallation reactions and cosmic ray data. 
Reactions of importance for the Li isotopes are protons and helium nuclei on C, N, and 
nuclei, and there is a contribution from helium-helium fusion. Considering recent results for 
the cosmic ray flux, composition, and propagation and incorporating the effects of ionization 
of interstellar gas, the computed ''Li/^Li ratios are 1.5-1.7 and 1.35 it 0.40 ^, values equal 
to the lower of the ratios found for o Per. Our conclusion is that the cloud material has 
been subjected to a severe flux of energetic particles. Such a flux has not affected ( Per's 
cloud. The line of sight toward o Per lies closer than that of ( Per to a site of star formation, 
IC 348, where massive stars could supply energetic particles via supernovae or stellar winds. 
Confirmation is provided by a study of the interstellar OH and HD chemistry of these sight 
lines that showed the cosmic ray fiux toward o Per was approximately an order of magnitude 
higher i'''^®. 

An earlier study ^° suggested a low ''Li/^Li ratio for a cloud toward ( Oph. Using the 
stronger K I line at 7699 A as a velocity template, isotope ratios of 8.8 it 0.8 it 1.4 and 
1.4to'5 ±0.6, where the second uncertainty represents systematic effects, were obtained. One 
concern with this two-component model for ( Oph is that the main velocity component is 
known to be 2 components separated by 1 km s~^ Another concern is that there is no 
evidence for molecular absorption from the second component toward ( Oph, unlike the 
situation in other directions Because this earlier study probably used an inappropriate 
template, the difference from the Solar System ratio of 12.3 is at best suggestive. Our 
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measurements, acquired at higher spectral resolution and analyzed with a more appropriate 
template, offer the first clear and unambiguous evidence for freshly synthesized products of 
cosmic ray spallation. 

Derivation of the total interstellar abundance is not trivial because knowledge of the 
abundance for Li II, the dominant ion, and mechanisms for depletion onto interstellar grains 
is necessary. An estimate for Li II abundance from data on Li I requires the electron density, 
which can be approximated by the amount of C II, the most abundant element providing 
electrons in neutral interstellar gas. The resulting total Li/H abundances appear in Table 
2. Since a Li isotope ratio of ~ 2 toward o Per suggests newly processed Li, the comparable 
Li/H abundances are unexpected. The two abundances, however, are very similar to recent 
interstellar Li results, using the analysis described in the legend to Table 2 for consistency: 
13.1 X 10-^° toward ( Per 11.0 x 10"^° toward p Oph ^ 20.6-24.4x 10"^° toward ( Oph ^^-^o. 
For further comparison, the Solar System value ^ is 20 X 10~^°. Since a weighted interstellar 
average of C II was utilized in determining toward o Per, there is additional uncertainty in 
the derived lithium abundance. There also remains the poorly known correction for depletion 
onto grains. 

Another observational measure may shed light on this inconsistency, the K/Li ratio. It 
is not dependent on electron density or the strength of the radiation field, and the amounts 
of depletion onto grains are similar The inferred K/Li ratios displayed in Table 2 for 
the clouds toward o Per and C, Per show no clear trend and are similar to the Solar System 
ratio ^ of 66 lb 8. Again, no indication for fresh Li is evident. 

Although the measurements are as yet few in number, a reasonable assumption based 
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on isotopic Li measurements for several clouds, local stars and the Solar System is that 
''Li/^Li ~ 10 once characterized o Per's clouds. If gas with such an isotopic ratio is exposed to 
energetic particles, fresh ''Li and ^Li are synthesized but little lithium is destroyed. Dilution 
of ''Li in this fashion would necessarily not attain the observed low isotopic ratio until the 
total lithium abundance had been raised about an order of magnitude. Neither the Li 
abundance nor the K/Li ratio in o Per's clouds are, however, unusual; Li appears to be 
less abundant. This suggests that the energetic particles operated on gas highly depleted in 
lithium: ejecta from massive stars - winds and supernovae - may be the playground for the 
energetic particles. 

Models ^'' of superbubbles, which contain ionized gas blown by supernovae in star forming 
regions, show that energetic particles synthesize about 10^^ Li atoms. At the end of the 
superbubble's existence, these Li atoms are diluted with swept up ambient material at which 
point the Li/H ratio is 100 times less than observed toward o Per. The clouds toward o Per, 
if spherical, contain 3 X lO"** Li atoms in 3 solar masses of gas. If the superbubble is itself 
roughly spherical, the total Li content inferred from the o Per sight line is close to the 
predicted 10^^ atoms. The superbubble scenario may account for our result provided that 
dilution of its products with ambient gas has not proceeded to completion. Sight lines 
containing partially mixed products from a superbubble like that toward o Per may be rare. 

The presence of deuterium as HD molecules toward o Per ^® implies that the gas within the 
superbubble cannot have been totally depleted in Li initially. The D-containing gas maybe 
largely in the cloud for which we find ''Li/^Li ~ 4. Synthesis of D within a superbubble is 
negligible. Existence of D may warrant consideration of alternative hypotheses, including 



6 



(i) isotopic fractionation of Li in diffuse clouds and (ii) the idea that pristine interstellar gas 
has ''Li/^Li ~ 2, the value for cosmic ray spallation, and varying degrees of contamination 
with ejecta from Li-rich red giants drive the ratio to the higher values seen in ( Per's clouds 
and elsewhere. 
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Table 1. Results of Profile Fitting. The synthesis included hyperfine structure for both 
^Li and ''Li The hyperfine structure splitting of K I A4044 was found to be negligible 
(~ 0.002 mA) and was not included in the analysis. The weaker member of the K I doublet is 
at 4047 A and thus is not part of the analysis. For oscillator strengths, we adopted tabulated 
results for each hyperfine transition in Li I and for the fine structure line of K I The 
profiles of the bluer ''Li I fine structure line and K I line were fitted first to obtain the 6- value 
and velocity of the line in the local standard of rest (u/sr), as well as column density (N) for 
each species. The 6- value, which is dominated by turbulent motion, and velocity were then 
fixed and used in the fit of the Li I profile. The uncertainties in N were determined from 
the uncertainties associated with the amount of absorption, while for the 6- value and u/^r, 
the results for the two fits provide a measure of their uncertainties. The final results toward 
each star are weighted averages determined by the separate fits. The results for component 
1 toward o Per are listed above those for component 2. The larger 6- value for the ''Li fit 
toward ( Per arises in part from the November data having a slightly larger instrumental 
width. For the weak lines examined here, the difference in 6- value produces a 10% effect in 
column density and in the inferred isotope ratio. 
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Parameter 



6- value (km s ^ ) 



visr (km s ^ ) 



NCU I) (10« cm-2) 



iV(6Li I) (10« cm-2) 



N{K I) (10^^ cm-2) 



o Per 

''Li lit K lit average 

2.1 2.1 2.1 

0.6 0.8 0.7 

6.5 6.5 6.5 

3.5 3.3 3.4 

26 ± 2 26 ± 2 26 ± 1 

5.6 ±0.6 5.3 ±0.8 5.5 ± 0.5 

7.0 ±1.5 7.4 ±1.5 7.2 ±1.1 

3.2 ±0.6 3.2 ±0.8 3.2 ± 0.5 

9.4 ±0.3 9.4 ±0.4 9.4 ± 0.2 

1.8 ±0.2 1.9 ±0.2 1.9 ±0.1 



C Per 

''Li lit K lit average 

1.9 1.4 1.7 

6.5 6.5 6.5 

37 ± 2 33 ± 2 35 ± 1 

3.7 ±1.4 3.0 ±1.1 3.3 ±0.9 

7.6 ±0.4 7.0 ±0.4 7.3 ±0.3 
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Table 2. Results for the '^Li/'^Li Ratio, the Total Li/H Abundance and the K/Li 
Abundance Ratio. For our estimates ol the total Li/H abundance, the electron density 
(ne) is obtained from N{C II), the total proton column density [iV(H) = iV(H I) + 2iV(H2)], 
and the gas density (n). Since no precise A^(C II) measurements exist for the direction toward 
o Per, the weighted mean interstellar ratio 2° of N{C II)/iV(H) = (1.42 ± 0.13) X 10"'* was 
utilized, yielding N{C II) = (2.2±0.8) x 10^^ cm'^. The value of N{C II) is (1.84±0.32) x 10^^ 
cm~2 toward ( Per '^^ . The columns iV(H) ^^'^"^ are (15.2 it 5.3) X 10^° cm"^ toward o Per 
and (15.8 lb 4.7) x 10^° cm"^ toward ( Per. Chemical models for carbon-bearing molecules '^'^ 
indicate values for n of 800 and 700 cm~"^ toward o and ( Per, respectively. The total Li/H 
abundance is derived from {[A^(''Li I)-|-A^(^Li I)]/A^(H)} X [G / (ane)], where a value of 41 
was used for [G/a)Li, the photoionization rate to recombination rate coefficient '^^ . The 
ionization corrections for the two clouds toward o Per are assumed to be the same. The 
K/Li ratios are based on column densities in Table 1 and a value of 9.4 for [G/a)K For 
the ''Li/^Li and K/Li ratios toward o Per, component 1 is the first entry and component 2 
is the second. 



Parameter 



o Per 



C Per 



Li/^Li ratio 



3.6 ±0.6 



10.6 ±2.9 



1.7 ±0.3 



Li/H abundance 



(9.8 ±3.5) X 10- 



(12.2 ±2.2) X 10- 



K/Li abundance ratio 



65 ±4 



42 ±2 



48 ±4 
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Figure Caption 



Fig. 1. Interstellar spectra of K I (a and c) and Li I (b and d) toward ( and o Per. The 
spectra were reduced with NOAO SUN/IRAF. After correcting for the bias and removing 
scattered light, the stellar images were flat fielded to remove pixel-to- pixel variations in 
sensitivity. Extracted spectra were produced by summing over pixels perpendicular to the 
dispersion direction and were placed on a wavelength scale with spectra from a Th-Ar hollow 
cathode taken throughout the night. Individual spectra were Doppler corrected and coadded, 
yielding a signal-to-noise ratio per pixel of approximately 2500:1 for each stellar spectrum 
based on the peak-to-peak variation in the continuum. The data points are filled circles. 
The dot-dashed line is the fit to the line profiles using the K I line for the velocity structure. 
The solid line gives the difference between data and fit; the offset is 1.002. The vertical lines 
in the K I plot for o Per show the positions of the two velocity components; a single velocity 
component is seen toward ( Per. The vertical dashes in the Li I plot for ( Per indicate 
positions of the ''Li and ^Li fine structure lines. In the absence of ^Li, an unsaturated 
''Li feature from a single cloud is a doublet with the red component half the strength of 
the blue component. A ^Li contribution augments the strength of the red component and 
contributes a weaker redder component. As long as the ''Li/^Li ratio is not small, the Li I 
feature is dominated by the 2:1 ratio of the blue to red component of (primarily) ''Li. For 
o Per, however, the red component approaches the strength of the blue component. Since 
K I A4044 shows that this enhancement cannot be due to a red-shifted ''Li component, the 
clouds toward oPer must be rich in ^Li. 
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